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In this article we present a concrete proposal for spin squeezing the ultracold ground state polar 
paramagnetic molecule OH, a system currently under hne control in the laboratory. In contrast 
to existing work, we consider a single, non-interacting molecule with angular momentum greater 
than 1/2. Starting from an experimentally relevant effective Hamiltonian, we identify a parameter 
regime where different combinations of static electric and magnetic fields can be used to realize the 
single-axis twisting Hamiltonian of Kitagawa and Ueda [M. Kitagawa and M. Ueda, Phys. Rev. 
A 47 , 5138 (1993)], the uniform field Hamiltonian proposed by Law et al. [C. K. Law, H. T Ng 
and P. T. Leung, Phys. Rev. A 63 , 055601 (2001)], and a model of field propagation in a Kerr 
medium considered by Agarwal and Puri (G. S. Agarwal and R. R. Puri, Phys. Rev. A 39 , 2969 
(1989)]. To support our conclusions, we provide analytical expressions as well as numerical calcula¬ 
tions, including optimization of field strengths and accounting for the effects of field misalignment. 
Our results have consequences for applications such as precision spectroscopy, techniques such as 
magnetometry, and stereochemical effects such as the orientation-to-alignment transition. 
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I. INTRODUCTION 

The ultracold OH molecule is a versatile platform for 
precision measurements. Its experimental appeal lies in 
the fact that its ground state is polar as well 

as paramagnetic, and is therefore readily manipulated in 
the laboratory with the use of electric and magnetic fields 
[IHi. Applications explored thus far include precision 
^ectroscopy in and quantum information processing 
Q, in addition to studies of cold chemistry [9I-[K^ and 
quantum degeneracy (l^ . 

In this article, we consider spin squeezing of the OH 
molecule. Our discussion occurs in the context of the 
Heisenberg uncertainty relation [T^ 

> l(J,)]/2, (1) 

between the three components of the angular momentum 
operator J. Spin squeezing refers to a situation where 
the fluctuation in one of the components, say AJ^,, is 
reduced to below the standard quantum limit ^J\{Jz)\/2. 
Of course, the fluctuations in Jy increase correspondingly, 
in order to maintain the relation of Eq. O- 

Spin squeezing constitutes a technique of interest at 
the frontiers of precision measurement and has appli¬ 
cations in spectroscopy, magnetometry, metrology, the 
detection of particle correlation and entan glem ent, and 
quantum computation and simulation (|I5l. Il6l| and ref¬ 
erences therein). The pioneering start to spin squeez- 
ing was provided by the work of Kitagawa and Ueda 
Hill and Wineland et al. , who considered the 

squeezing of collective atomic spins, and was followed by 
many investigations ( (2ll - l59jl . for example). We empha¬ 
size that all of this work relates to ensembles of correlated 


(pseudo)spins, represented by atoms (in a thermal vapor 
or a degenerate gas), nuclei [60j , or molecule s 16111 , where 
the collective spin can assume a high value [62j| . 

More recently, squeezing has also been considered for 
single (atomic or nuclear) spins, or equivalently, for an 
uncorrelated ensemble of suc h sp ins. Experiments have 
been carried out, using spin 3 and spin 7/2 states 
of Cesium atoms. Possibilities also exist for Dysprosium 
which offers states with spins from 8 upto 12.5. The¬ 
oretical calculations of single-spin squeezing have been 
presented as well [b^. The squeezing in this case is lim¬ 
ited by the much smaller angular momentum available. 
Nonetheless, some effects can be experimentally relevant, 
for example for precision spectroscopy [l^. Motivated 
by such a consideration, we consider spin squeezing of a 
single OH molecule. Interesting applications to the mea¬ 
surement of magnetic fields also seem possible [63 - [70l| . 
especially following recent discussions of magnetom etry 
using single SiC spins with angular momentum 3/2 [Tlj; 
ground state ultracold OH also carries rotational angular 
momentum 3/2, and is sensitive to magnetic fields via the 
Zeeman shift. As well, spin squeezing is related to the 
alignment-to-orientation transition |65l ]. and we expect 
this perspective will be of relevance to our work, given 
the recent interest in the stereochemical properties of the 
OH molecule [t^. We emphasize that compared to pre¬ 
vious proposals for spin squeezing molecules (^ . which 
employ spin 1/2 interacting molecules squeezed collec¬ 
tively, we consider a single noninteracting molecule with 
angular momentum 3/2. 

In present-day laboratories, ultracold OH molecules 
are typ ically confined in magnetoelectrostatic traps ii 
iTa uM ■ We therefore consider spin squeezing enabled 
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FIG. 1: (Color online) Schematic of the diatomic OH molecule 
in electric E and magnetic B fields crossed at an angle 9. 


by these readily available static electric and magnetic 
fields. The advantage of using static, rather than opti¬ 
cal or microwave [l5j, fields for squeezing is that damp¬ 
ing and decoherence due to spontaneous emission can be 
avoided. Our starting point will be an effective eight¬ 
dimensional matrix Hamiltonian that has been shown 
to model recent OH experiments quite well , 

and has also been diagonalized analytically . We 

demonstrate that for Zeeman and Stark shifts small com¬ 
pared to the Lambda-doublet splitting of the OH ground 
state, this Hamiltonian can yield spin-squeezing of the 
types considered by Kitagawa and Ueda [l^, Law et al. 
[ 3 ^ and Agarwal and Puri [2lj. We provide analytic and 
numerical results, discuss the optimization of the field 
strengths, and include the effects of field misalignment 
in our treatment. 

The remainder of this article is arranged as follows. 
Section El presents the derivation of the spin squeezing 
Hamiltonian, Section uni discusses the dynamics of the 
squeezing parameter. Section HVl addresses the detection 
of the proposed squeezing, and Section |V] supplies a con¬ 
clusion. 


volving only eight quantum states lzl,[ 2 i- The domain of 
validity of this Hamiltonian and the details of the states 
involved can be found in several articles I, El 111, [ 73 , 
and will not be repeated here. This effective Hamiltonian 
was recently diagonalized analytically following the de¬ 
tection of an underlying chiral symmetry (75l | . In the pro¬ 
cess of identifying that symmetry, it was found that the 
OH matrix Hamiltonian could be re-expressed in terms of 
two interacting spins (to avoid notational clutter we set 
= 1; to make contact with laboratory parameters, req¬ 
uisite factors of h can be supplied to any formula in this 
article by inspection, see the example supplied below), 

Hm = — Acr^ — BJz + E(Jx {Jz cos 9 — Jx sin 9 ), (2) 


where the constants are given by 



4^3 H 
5 


.E = 


2fJ,eE 

5 


( 3 ) 


A being the Lambda-doublet splitting, B and E the uni¬ 
form magnetic and electric fields, respectively, /is the 
Bohr magneton and /ie the electric dipole moment of the 
OH molecule. In Eq. m,e is the angle between the elec¬ 
tric and magnetic fields (see Fig. [T|), and a = 1/2 and 
J = 3/2 are the two interacting spins. The spin cr is 
a pseudo-spin, with the spin projections cr^ = — 1 and 
CTz = 1 corresponding to the two Lambda-doublet mani¬ 
folds of parity e and /, respectively. On the other hand, 
J corresponds to the rotation of the molecular axis, and 
the values Jz = ±1/2, ±3/2 indicate the various projec¬ 
tions of J on the laboratory z axis, which is chosen for 
convenience and without loss of generality to be along 
the magnetic field. The matrix representation of EIm is 
reproduced in the appendix for the reader’s convenience, 
see Eq. ((151) . 

At this point the insertion of factors of h in Eq. @ 
may be considered as a useful exercise. We begin with 
the observation that EIm has units of energy. Since the 
Pauli matrix ax is dimensionless, A needs to have units 
of energy, and thus must be multiplied by h. On the 
other hand, has units of h. Therefore B needs to have 
units of frequency, and thus the right hand side of the 
corresponding expression in Eq. ([5]) must be divided by 
an h. Other formulas in this article can be handled in a 
similar manner. 


II. DERIVATION OF THE SPIN SQUEEZING 

HAMILTONIAN The four dimensional squeezing Hamiltonian 

A. The eight dimensional effective Hamiltonian generation of spin squeezing requires the presence 

of a nonlinearity, represented to lowest order by a term 
quadratic in one of the spin operators in the relevant 
Hamiltonian E3l- However, Eq. m is linear in each of 
the spin operators. It may not be readily obvious that the 


Several experiments on ultracold OH molecules in 
crossed electric and magnetic fields, see Fig.I, have been 
successfully modeled using an effective Hamiltonian in- 
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coupled dynamics of the two spins can lead to squeezing, 
and it would assist our intuition if Hm could be reduced, 
even if in some restricted regime, to the form of a spin 
squeezing Hamiltonian familiar from the literature. In or¬ 
der to effect such a reduction, we proceed by identifying 
a regime in which the pseudospin 1/2 can be eliminated 
adiabatically. As we show below, this results in an effec¬ 
tive spin-squeezing Hamiltonian for the spin 3/2 degree 
of freedom. A similar procedure has been used earlier in 
the case of nuclear spin squeezing; in that case however 
the eliminated spin is real (^ . The derivation begins 
with the Heisenberg equations implied by Eq. ([2]) for the 
system variables 

= 2Aay, ( 4 ) 

&y = —2Aax — 2Eaz {Jz cos 9 — sin6>), (5) 

az = 2Eay {Jz cos 9 — sin0), (6) 

Jx = BJy — Ecos9axJy, (7) 

Jy = —BJx+ E(Tx {Jx cos 9 + Jz sin 9) , (8) 

Jz = —Esm9axJy ( 9 ) 

We now consider the regime 

A>E,B, (10) 


i.e. such that the Lambda-doublet splitting is larger than 
the Stark as well as Zeeman shifts. In our calculations be¬ 
low we will work with A/A ~ 0.25, which corresponds to 
an electric field of 100 V/cm, the lower bound stipulated 
by the value of stray fields currently affecting experiments 
(74l| . We will also consider B/A < 0.1, corresponding to 
a magnetic field of 20 G, which can easily be achieved 
exp er iment ally. 

In the regime indicated by Eq. m, Gx and Gy will 
vary at the high rate ^ 2A and may be thought of as 
the ‘fast’ variables of our problem. We will see below 
that with this assumption Gz ceases to be a dynamical 
variable, i.e. takes on a constant value. On the other 
hand, Jx,Jy and Jz^ varying at the low rates ~ B and 
^ E (to first order), may be thought of as the ‘slow’ 
variables of the problem. We now adiabatically eliminate 
the fast spin variables Gx and Gy in Eqs. 0-® by setting 
the time derivatives in those equations equal to zero. The 
resulting solutions are 

EC 

Gx = {JzCOs9 - JxS\n9), ( 11 ) 

a, = 0, (12) 

a. = C, (13) 


where we have retained C for generality. Note that the 
adiabatic solution Gy = Q forces Gz to be a constant [see 


Eq. (0], as mentioned above. We now use the solutions 
of Eqs. (III1)-(II31) to adiabatically eliminate the fast spin 
1/2 degrees of freedom from Eq. ([2]). Dropping a con¬ 
stant term proportional to A, we obtain the adiabatic 
Hamiltonian for the spin 3/2 variables 

ce‘^ 

Ela = -BJz -^ {Jz COS 9 - Jx sin 0)^ , (14) 

which clearly contains terms nonlinear in the angular mo¬ 
mentum components, and can therefore enforce squeez¬ 
ing. 

Below, we will show that some standard squeezing 
Hamiltonians can be recovered from Eq. o, and also 
explore the squeezing effects of this Hamiltonian for ar¬ 
bitrary B, E and 9, within the limits prescribed by the 
validity of the effective Hamiltonian [Eq. ([2])]. Subse¬ 
quently, we will compare these results with numerical 
calculations based on the full eight-dimensional Hamilto¬ 
nian of Eq. which accounts for non-adiabatic effects 
and thus allows us to examine the validity of the adi¬ 
abatic elimination process used in arriving at Eq. m- 
We note that Eq. does not stipulate any relation¬ 
ship between E and B, other than that they both have 
to be smaller than A. We will use this flexibility below 
to adjust the field magnitudes to optimize squeezing. 


III. SPIN SQUEEZING DYNAMICS 

A. Spin squeezing using the four-dimensional 
adiabatic Hamiltonian Ha 

In this section we consider spin squeezing using the 
four-dimensional adiabatic Hamiltonian of Eq. (fHl) . We 
provide analytic and numerical results, as appropriate. 


1. One-axis twisting 

For B = 0,9 = OjC = 1, Eq. (fHl) gives the Kitagawa- 
Ueda Hamiltonian [Iq . 


Hku = kJ'^, 


where 



(15) 


(16) 


The theoretical analysis for spin squeezing using 
was first provided by Kitagawa and Ueda [l^ , and more 
recently by Rochester et al.[6^. Nonetheless, we restate 
the procedure here, in order to compare with fully nu¬ 
merical calculations to be presented later. The recipe for 
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the analysis goes as follows. Using the four-dimensional 
matrix for the J = 'i 12 operator we write the ma¬ 

trix form for Hku- We then obtain the time evolution 
operator 

Uku = (17) 


which is also a four-dimensional matrix. We can then 
obtain the time evolution of any observable O by using 
the relation 0{t) = U]^]jOU ku■ In this way we find the 
matrix forms of the observables Jx{t) and Jy{t). In order 
to find the most suitable axis for spin squeezing, it is 
convenient to consider a further rotation by an angle n 
about the x axis [l^, i.e. 

(18) 

and so on for other observables. The operators Jx{t) and 
Jx{t) are unaffected by this rotation, of course. 

For our initial state, we choose, following Kitagawa and 
Ueda, the (coherent) stretched state along the x axis [l8l| . 


NW = |j = 5.M = i), = 2-»/=y: 


i! 

1 2’ 2 ’ 

k=0 ^ ' 

( 19 ) 

which we have expanded in the 2 basis on the right hand 
side. We note that the current level of control over the 
OH ground state manifold should readily allow this state 
to be prepared M- Using the initial state \i}xu 
find the expectation values {Jy,n{t)) , {Jy .^{t)) , etc., and 
thence the variances such as 

(AJ,.„(t))^ = (4„(t))-(J,,„(t))^. (20) 

This procedure yields the relevant quantities @ 
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{Ut)) 

(AJ,,„(f))2 

where 


■ cos^ kt. 


( 21 ) 


^ M VM2 + iV2 , 

1 -f -f-^-cos { 2 ti -\- 2J) 


, M y/JPTW 
1 -I- -^- COS (2n -I- 25) 


( 22 ) 


M = 1 — cos2kt, N = 2sm2kt, 5 = - tan ^ 

2 


(23) 

As pointed out earlier the y quadrature is max¬ 

imally squeezed for cos (2nopt + 2^) = — 1, i.e. along the 
axis 


^opt — 2 


(24) 


in the y — z plane. The amount of squeezing is quantified 
by the parameter introduced by Wineland et al. [2^ 




y,n 




and, analogously, 


\{um 


(25) 


(26) 


Squeezing occurs when either or ^z,n is less than one. 
We have shown a plot of ^y,riopt in Fig- 111 The blue (solid) 
line is the analytical prediction of Eq. (1^51) and shows that 
the y quadrature is squeezed periodically. As found by 
Rochester et al. [^, for J = 3/2 the minimum value 
of the squeezing parameter is — 0.75, consistent 

with Fig. [21 The red (dashed) curve in Fig. [2] corresponds 



kt 


FIG. 2: (Color online) Plot of a function of the 

dimensionless time kt where k is dehned in Eq. (1161) . The pa¬ 
rameters are A = 1.66 GHz, E = 100 V/cm and 5 = 0, im¬ 
plying K = 48 MHz. Squeezing occurs when ^y,nopt < 1- The 
blue (solid) curve is the analytical prediction from Eqs. (1241) 
and (12511 based on the adiabatic approximation of Eq. as, 
with E/A ~ 0.25. The red (dashed) curve is the numerical 
calculation using the eight dimensional effective Hamiltonian 
of Eq. (|2|. As can be seen, nonadiabatic effects significantly 
change the magnitude and periodicity of the squeezing. 

to a numerical calculation of squeezing using the full ef¬ 
fective Hamiltonian of Eq. ©, which incorporates the 
effects of nonadiabaticity, see Section IHI HI below. As 
can be seen, nonadiabatic effects change the magnitude 
and periodicity of squeezing quite significantly. This may 
be expected, as our adiabaticity parameter E/A ^ 0.25 
is quite large. A plot of ^z,nopt has been shown in Fig. O 
The analytical prediction of Eq. (E51) , represented by the 
blue (solid) curve, implies that the shown z quadrature 
is never squeezed. Interestingly, the inclusion of nonadia¬ 
batic effects, shown by the red (dashed) curve, introduces 
some squeezing of this quadrature. 
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FIG. 3: (Color online) Plot of ^z.nopt as a function of the 
dimensionless time kt where k is defined in Eq. (I16II . The pa¬ 
rameters are A = 1.66 GHz, E = 100 V/cm and 5 = 0, im¬ 
plying ft = 48 MHz. Squeezing occurs when ^z,nopt < 1- The 
blue (solid) curve is the analytical prediction from Eqs. (1241) 
and (1261) based on the adiabatic approximation of Eq. 113, 
with E/A ~ 0.25. The red (dashed) curve is the numeri¬ 
cal calculation using the eight dimensional effective Hamilto¬ 
nian of Eq. As can be seen, nonadiabatic effects actually 
introduce some squeezing into the otherwise ‘anti-squeezed’ 
quadrature. 


2. Uniform field Hamiltonian 
In this section we consider spin squeezing in the pres- 


ence of a magnetic field. 

For B ^ 0, 

and 6 = 7r/2, 

Eq. ( 1151 ) takes the form 



Hlnl = 

— BJz ^Jx’ 

(27) 


This Hamiltonian is, to within a unitary anti-clockwise 
rotation of n/2 around the y axis, the same as proposed 
earlier by Law, Ng and Leung [^. [78l - [80l| . These authors 
suggested adding the uniform field term to the Kitagawa- 
Ueda Hamiltonian as it yielded greater squeezing for 
longer times. Analytic solutions to the spin squeezing dy¬ 
namics of Eq. (1271) are not available for arbitrary spin J. 
Generally, analytic results can be found only for J < 2, 
since the eigenvalues of the Hamiltonian are required for 
the calculation, which can be found in closed form only 
for matrices of dimension 5 or lower. Some results for 
the case where J = 3/2 is a collective spin have been 
published in the literature Hill. We provide addi¬ 
tional expressions in order to discuss the details of our 
problem. 

To determine the squeezing, we follow a procedure sim¬ 
ilar to that of Section IHI A H but consider the initial 
stretched state along the z direction 


and determine the squeezing about the x and y axes. We 
find 


(Mt)) 

{{Ajyit))y 


{Jy{t))=0, 



(29) 

(30) 

(31) 
,(32) 


where 

P = \/B‘^ - Bk + k?. (33) 


Using Eqs. (|^ - (IM1) . we can find the squeezing parame¬ 
ters 


and 






(34) 




/- ((AJy(t))) 


(35) 


Plots of the squeezing dynamics are shown in Fig. HKa) 
for the same parameters as in Fig. ([2]) but with B = 20 
G. The dashed curves represent the analytical results for 
[Eqs. (IMl) ] (blue) and ^y [Eq. (1551) ] (red). The solid 
curves are the corresponding numerical calculations us¬ 
ing Hm from Eq. ([5]). It can be seen from the numer¬ 
ically calculated red (solid) curves in Fig. |4][a) that the 
Law-Ng-Leung approach is quite robust to nonadiabatic 
effects. In Fig. |4| we have chosen a magnetic field opti¬ 
mized using the following procedure. From the analytic 
result of Eq. (1351) , it can be seen that at multiples of the 
time ts 


ts = 


TT 

Ip’ 


(36) 


the squeezing parameter for the y quadrature attains an 
extremum value given by 


= ts) 


2^{r^-r + l) (r2-2r-b2) 
2 r2 -2r+l 


where the dimensionless ratio 


(37) 


r = 



(38) 


The denominator of ^y(t = ts) vanishes only for the com¬ 
plex values r = (l±i)/2 , which are excluded by experi¬ 
ment. Thus, stays finite as r varies, as can be seen 


I*) LNL ~ 1) 


(28) 
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FIG. 4: (Color online) Plots of as functions of the dimensionless time Pt where P is defined in Eq. (I33|l . The parameters 
are A = 1.66 GHz, E = 100 V/cm, H = 20 G, implying P = 144 MHz and (a) 6 = 90° (b) 6 = 85° (c) 6 = 80° (d) 61 = 75°. 
Squeezing occurs when < 1- The dashed curves are the analytical predictions for [Eq. (1341) ] (blue) and ^y [Eq. (1351) 1 (red). 
The solid curves represent numerical calculations carried out using the eight dimensional effective Hamiltonian of Eq. ([2|. As 
can be seen by comparing with Fig. O the squeezing is quite robust to nonadiabatic effects. Also, deviations of about 5° away 
from the nominal value of 90° for field alignment do not affect the squeezing greatly. 


from Fig. [SJ Differentiation of Eq. (1571) readily yields 
a minimum value = ts) — 0.8, which occurs at 

r ~ 3.3, indicated by the vertical line in Fig. [5] This 
optimized value of r corresponds to the magnetic field 20 
G used in Fig. IHa). For simplicity, unlike in Ref. [s^, 
we have not found the time-dependent axis of optimum 
squeezing, which can lead to even better squeezing than 
we have presented. 


3. General case 

In this section we consider spin squeezing as a function 
of the angle between the electric and magnetic fields. The 
interest in this degree of freedom arises from the neces¬ 
sity of accounting for possible misalignments between the 
electric and magnetic fields in the laboratory, e.g. away 
from the nominal value oi 9 = 7r/2 for i/iATL, see Sec¬ 
tion [TITXJl While this case can be solved analytically as 
well, the expressions are lengthy and we provide numer¬ 


ical solutions instead. The Hamiltonian is 

Hg = —BJz + k (Jz cos 9 — Jx sin0)^ . (39) 

Before we proceed further, we mention that this Hamil¬ 
tonian, when rotated anti-clockwise by an angle 9 about 
the y axis yields 

Hg = = —B (Jz cos 9 — sin0) -I- kJ^, 

(40) 

which is of the form considered earlier by Agarwal and 
Puri for arbitrary J [2l[. 

Numerical plots are presented in Fig. |4]for (b) 9 = 85° 
(c) 9 = 80° and (d) 9 = 75°, respectively. The dashed 
curves are the numerical implications of Eq. (13911 for 
[Eq. (IMl) ] (blue) and [Eq. (1551) ] (red). The solid curves 
are the corresponding numerical calculations starting 
from Eq. As can be seen, the squeezing is quite 

robust to field misalignment. Only at about 9 = 75° 
does the pattern change noticeably from that at 0 = 90°. 
Interestingly, while misalignment degrades squeezing in 
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From here the expectation values of the relevant opera¬ 
tors can be calculated, such as 

(J,(t))=Tr3/2[J./93/2(i)], (45) 


and therefore also the squeezing parameters, as for ex¬ 
ample in Eq. (1341) . It may be useful to mention that for 
Section IIII A II the initial state is now written as 

I*)xc7 = ^|0,0,0,0,1,V3,V3,1), (46) 


and for Sections [ 


r 

FIG. 5: (Color online) Plot of ^^(t = ts) [Eq. (I^ j versus 
r [Eq. (I38l) [. The vertical line corresponds to r ~ 3.3, at 
which ^y(t = ts) has a minimum, corresponding to maximum 
squeezing. 

the y quadrature, it does not correspondingly introduce 
squeezing in the x quadrature. 

B. Spin squeezing using the eight-dimensional 
effective Hamiltonian Hm 


l*)LArL = 10,0,0,0,0,0,0,1). (47) 

Comparison of the four and eight-dimensional results can 
be seen in Figs.lJHl and has been described in the preced¬ 
ing sections. These plots imply that while non-adiabatic 
effects can be important, the adiabatic approximation 
captures the basic spin dynamics reasonably well, and 
the intuition afforded by this approximation about the 
existence of squeezing is well founded. Better agreement 
can be produced by assuming a smaller adiabaticity pa¬ 
rameter, but currently this is limited by stray electric 
fields in the experiments. The key point is that even in 
the presence of non-adiabaticity, squeezing exists. 


In this section we describe our spin squeezing calcu¬ 
lations using the full eight-dimensional Hamiltonian Hm 
of Eq. (I2|). This enables us to quantitatively examine 
the adiabatic approximation made in deriving the four¬ 
dimensional spin squeezing Hamiltonian Ha of Eq. da. 
The calculation proceeds as follows. We begin with 
the eight-dimensional matrix representation of Hm from 
Eq. (H51) . Although the treatment can be carried out an¬ 
alytically, the expressions are very long, and we calculate 
instead numerically the time evolution operator 

Um = (41) 

which is also an eight-dimensional matrix. Starting 
from the initial state |'!/’(0)), represented by an eight¬ 
dimensional column matrix, the state vector \ip{t)) at 
any later time t is 


m)) = Um \m) • (42) 

The density matrix of the full system can then be easily 
foun 


p(f) = IV’(t)) (V’(i)l ■ (43) 

The reduced density matrix of the spin 3/2 system can 
then be found by tracing over the spin 1/2 degrees of 
freedom, i.e. 


P3/2(i) = Tri/2 [p(f)]. (44) 


IV. DETECTION OF SPIN SQUEEZING 


The proposed spin squeezing may be detected by per¬ 
forming quantum tomography on the OH ground state 
manifold, which would yield the density matrix, from 
which squeezing information can readily be extracted. 
Such procedures can be carried out for OH in analogy to 
experiments performed earlier on atomic and nuclear 
64| systems. In the laboratory, the predicted squeezing 
will be degraded by damping and noise, due to molec¬ 
ular collisions and trap loss. In the present work, we 
have justifiably neglected these effect^ since they occur 
at typical rates of Hz (collisions) [1, 0, 0 or KHz (trap 
loss) [Z^, while squeezing is generated at frequencies of 
MHz [Fig. [5]. 


V. CONCLUSION 

We have proposed a scheme for spin squeezing the ul¬ 
tracold OH molecule, in the context of ongoing experi¬ 
ments. Production of such nonclassical states is expected 
to be useful for spectroscopy, magnetometry and stere¬ 
ochemistry. We have identified an accessible parameter 
regime where single-axis twisting as well as uniform field 
squeezing can be implemented with the use of static fields 
only. Since we do not propose to use optical or microwave 
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fields, our scheme is free from damping and decoherence 
due to spontaneous emission and optical pumping. In 
our analysis, we have shown how to optimize the field 
values and also investigated the effect of field misalign¬ 
ment on the squeezing. To our knowledge, our work is the 
first concrete proposal for spin squeezing single noninter¬ 
acting molecules with angular momentum greater than 
1/2. We note that with the use of additional electric 
and magnetic fields, other spin squeezing Hamiltonians 
may also be realized using OH, such as the one proposed 
by Raghavan et al. [3l|. Also, a more accurate descrip¬ 
tion of the OH ground state can be reached by including 
more states in the Hamiltonian, accounting for fine and 
hyperfine structure and electric quadrupole interactions 


Sil- It would be interesting to investigate the effect 
of these additions on our results. Finally, our scheme 
can also be extended to the ground states of other polar 
paramagnetic molecules such as ^ 113/2 LiO and ^$2 CeO. 

We would like to thank S. Marin for useful discussions. 


VI. APPENDIX 

In this Appendix we provide, for the reader’s conve¬ 
nience the matrix representation of Hm from Eq. @ 

[zHil, 



0 

0 

0 

^fleE COS 9 

— ^HeE sin 9 

0 

0 \ 

0 


0 

0 

— ^^eE sin 9 

^fieE cos 9 

— j^eE sinS 

0 

0 

0 


0 

0 

— l^eAsin^ 

— ^HeE cos9 

— ^fieE sin 9 

0 

0 

0 


0 

0 

— ^fieEsm9 

— l/TeA cos 9 

l/TeU cos 9 


0 

0 

¥ - IfsB 

0 

0 

0 

— ^fieEsinO 

i/TeU cos 9 

— ^HeE sin0 

0 

0 

2 , , R 
— - g/isH 

0 

0 

0 

— |/ieAsin0 

— ^HeE cos9 
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